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Q^ , Abstract: The consideration of the N-body gravitational problem equations 

0^ ' can give to us some class of boundary- value problems defined on the "beem's" 

construction. One can considere it as weak or so-called finite element method's 

r^ ' approximation with the linear function's basis. So there are time and coor- 

5^ , dinate dependent beems as distances between points and unknown coordinate 

functions which are defined on the beems. The possible sample from this class 

OO , is simple Lalplas equation where potential energy are maximal. Some com- 

puter simulations can give to us symmetries on spherical coordinates (similar 
for well known mathematical structures - configurations like Pascal or Desargus 

f^ ' configurations). 

^SJ , But in other hand, the consideration of the natural mass distribution can 

give to us evidence of the physical reality of the similar structures on the star 

^^ , sphere. Moreover, there are coupled invariancy for geometry and arithmetics. 
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Introduction. 

So, it is well known that there is a structure of young O, B and AQ stars. 

Herschel (1847) was the first who noted and statistically convincing arguments 

were given by Guld (1879), that there is so called Guld zone. There was dis- 

H \ cussion of physical validity and stability- problem of so called O, B associations 

and superassociations [1], [2], [3]. At the last time appeared the results on the 
periodicity for the clusters of the galaxies [4] , [5] . Possibly there is mathemat- 
ical unity for some astrophenomena or in other words automodelity for mass 
distribution so J.Einasto wrote: ' If this refiects the distribution of all matter 
(luminous and dark), then there must exists some hitherto unknown process 
that produces regular structure on large scales' [5]. But word 'unknown' does 
not mean existence of really new physical law which take place only for large 
scales. We hope that there are mathematical methods yet and order for dis- 
tances smaller than cosmological sizes. We start our investigation from famous 
and old N-body problem [6], then for assemblies of bodies with low velocities 
we can find coherent structures. 



But there is similar visual coherentness, moreover automodel symmetry for 
some scales on natural coordinate systems ( spherical systems of coordinates 
related with momentum). 

N-body equations for x,y, z — 0. 

Let us consider N-body equations for N points with masses Mj on OXYZ 
coordinate system [7]: 

j=N' 

iij = ( ^ Mj{uj ~ Ui)/rfj), Ui = Xj, j/i, Zi (1) 

where r^ is distance between points i and j , sign ' noted so there is no case 
After decomposition we can write: 

3=N' o=N' 

u, = -u,iJ2 M, /r% ) + ( 5] uj Mj lr% ),u, = x,,y„Zi (2) 

Note, that right parts of (2) are equal to product of the positive matrix 
A{x,y^z) (by weak diagonal predominance) on the vector w = (xi,X2, --xn, 
yi,..yN,-.-.ZN)- 

Then let us consider (2) as finite approximation of some boundary-value 
problem on the "beem's" mesh where lenghtess of beems or in other words 
distances between points must be solution depended. 

We'll try find form from our finite-difference equations or 'revcrs motion' on 
well known transformations: 

Pt + ip4>x)x + ip(f>v)v = 0, (/)„ = (3) 

v{x, y){pt + {p(t>x)x + {P4>y)y) = 0, (4) 

{vpt - p4>xVx - p4>yVy)duj = (5) 

in 

For right part we can find function f{l) from: 

RightPart = / f{l)uividl (6) 

where L is beem in the distinguished moment of the time and integration are 
going along directions between points, v{l) is test function. 

Then we can take finite basis of functions Vi{l) (Figure 1) which related with 
points and which have finite region where v{l) is not equal to zero (so-called 



finite elements [8]). In (3) wc have: 



k=N 



i(o = Y. "'^^fc(^) (^) 



k=l 

After that we have that /(/) — Mj/ R^~ constant on every beem where R — rij. 
For left part we have: 

LeftPart - 2/K{ f v{l)uu/R - R/Guumdl) (8) 

where K is the number of the nodes which significant connected with the given 
one. 

In differential form and for the equal masses M we have: 

Mutt/R + {MR/6utu)i = {MyR^ui)i (9) 

where we take K ^ 2 without less of generality with respect to physical case. 

We can consider the row of equations where one of them is (6) and others 
are: 

M/Rutt = {M^/R^ui)i (10) 

{M/Rui)i = (11) 

where R is equal to small number when we take by unit size of our constel- 
lation and gravitational interaction between points defined on finite distance 
which are less than MR for some natural M . 

Let us consider (11). We can see so there are spherical symmetry for form 
(11) and hence there are some hope for symmetries on spherical representations 
on 'natural' ( that it to say sperical) coordinates. 

There are easy way for computer simulation. Let us consider volume G and 
points with equal masses where some points (set A) are into the volume G and 
other points (set B) are outside the volume G. The every point is surrounded by 
an P envelope of some (set P is not equal A+B) points. Note, near G boundary 
Pz have points from A and B too. 

Let's define iterative motion for every point z from volume G into mass 
center of P^. Then after some steps we will have solution for (11). 

We can see some examples of simulation on Figures 2,3,4 where P defined 
by choosen distance. 

By the way — it is simple and fast algorithm for Steiner problem [10] for 
two or three dimensional cases — Figure 5. 

Application 



Well known so equation (11) seems as fully unacceptable for real masses and 
formal solutions of (11) seems as very inconstant by time. 

But let us consider nearest and well known light galaxies — M31, M81, M82, 
M83, MlOl, NGC 5128, Dwingelool (data from nedwww.ipac.caltech.edu). 

We can see pictures — Figures 6,7,8,9 which are similar for our computer 
simulation. 

At next let us consider "histogram" for natural stars as distribution on the 
star's sphere or in other words as the distribution on the one of the planes:(^, b), 
{a,S) or (A,/3) on the one or two coupled rectangles of size 2tt * tt { there are 
periodicity). Also there is physical value for star which can be correlated with 
possible order. It is so-called visual brightness or in the order words order which 
was related with some degree of the M"^ /BP'^ where R is distance between star 
and view point [11]. Below are given the enumeration for 50 lightest stars [11]: 

1 aCma 2 aCar 3 aBool 4 aLyr 5 aCen 6 aAur 7 f30ri 8 aCmi 9 aOri 10 
aEri ll/SCen 12 aAql 13 aCru 14 aTau 15aSco 16 aVir 17 PGem 18 aPsa 
19 (3Cru 20 aCyg 21 aLeo 22 eCma 23 aGem 24 XSco 25 jOri 26 jCru 27 
PTau 28 pCar 29 eOri 30 aGru 31 eUma 32 (Ori 33 aUma 34 aPer 35 jVel 
36 eSgr 37 SCMa 38 rjUma 39 eCar 40 OSco 41 l3Aur 42 -/Gem 43 aTra 44 
6Vel 45 aPav 46 l3GMa 47 aHya 48 oGet 49 aAri 50 aUMi 

As an example, let us consider the stars with numbers 5 aGen, 6 aAur, 8 
aGmi, 12 aAql, 14 aTau, 15 aSco on (A, /3) coordinate plane ( or on cilinder) 
- Figure 10. 

It is an approximate example of graph of a new kind, similar for graceful 
graphs [12] and harmonious graphs [13], where arithmetic, geometry and topol- 
ogy are tightly coupled. 

We can give definition: 

each point i belongs to a straight line drawn through other two points j,k from 
the list, the sum of j and k being equal to 20. 
where is exception to the rule for point 5. 

Also we have on the {l,b) and {a, 5 ) planes. Figures 11, 12 — 

for every points i,k,j which belong to a one of the straight lines there is other 
straight line with points i, j, k, or i, j, 20-k or i,j-20,k-20 (on (l,b)) . 

Note that these are really the nearest extremally massive objects ( with 
exception of point 15 ). Moreover, there are the correlations with velocity vectors 



for presented configurations (or quasi-configurations) [14] in the all planes by 
time. 

Also we can construct similar figures for other viewpoints- view from aCmi 
give to us Figure 13, view from aCen Figure 14, view from aAql Figure 15. We 
can see so there are same numbers — 5, 6, 8, 12, 14, 15. 

Let us consider the stars with the numbers 5fc and stars nearest by same 
order. We have paradoxical meshes of the parallel lines for these stars- Fig. 
16 (A,/3) plane. Figure 17 {a,S ) plane, and Figure 18- for {l,b) planes (with 
parallelism on projective sense). Note that these are really extremally massive 
objects- more than 5 Sol's mass. 

With rare exception, most part of O, B and AQ stars (brighter than 6.0™) 
is concentrated near the straight lines drawn through pairs of the stars with 
numbers 2, 3, 5, 6, 7, 10, 11, 15- Figure 19; that is true for {l,b), {a, 6), (A, /3) 
coordinate planes. 

There is a substructure of stars of j3Per (EA) and (3Lyr (EB) types, which 
concentrate only to some of the lines in all the coordinate systems. It is an 
independent test, the brightness of these stars is 8 — 9™ on the average. K,M 
giants also gravitate towards these lines, but their positions are more smeared 
than those of O, B stars. 

Conclusion 

Perhaps phenomenon does not completely visual and there are order on 
native coordinate representation for equations (1) and (11). Maybe gravitational 
interactions for far and large masses are masked by heat motions of small stars 
and interstellar media. 
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Fig.l Example of nodes related linear function 
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Fig. 2. Simulation, 4 points, spherical plane, one on the viewpoint. 



100 




c 



c 



Fig. 3. Simulations result, 5 points, one on the viewpoint. 




10 




o • 



-Q 



O 



Fig. 5 Circles - 'boundary' points, filled circles-points 
in the 'volume', 2 and 3 D Steiner configurations. 
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Fig. 6. Light and nearest galaxies, lighter than 8.5™ -(l,b) representation. 
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Fig. 7. Light and nearest galaxies, lighter than 8.5"*, 
(A, /3) representation. 
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Fig. 8. (a, 6)- representation. 
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Fig. 9. Supergalactic plane representation. 
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Fig 10. 5- aCen, 6-aAur, 
8- alphaCmi, 12 -aAql, 14- aTau, 15- aSco (A,/3) plane. 



16 




17 




18 




Fig 13. Viewpoint is aCmi, 5-aBoo, 
6- PGem, S-aTau, 12- aCru, 14- aGem, 15- aSco 
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Fig. 14 Viewpoint is aCen, 5-aLyrfi- aAur,8-aOri, 
12- aCru, 14- aSco, 15- aVir 
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Fig. 15. Viewpoint aAql 
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Fig 16 (A, 13) "multiplied by 5" 
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Fig 17 (a, (* 


') "multiplied by 5". 
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3 Fig. 18 For first 50 stars on (l,b) 
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Fig 19 For first light 500 stars on a,/?. 
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